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Abstract. Set membership of points in the plane can be visualized by connecting 
corresponding points via graphical features, like paths, trees, polygons, ellipses. 
In this paper we study the bus embeddability problem (BEP): given a set of col¬ 
ored points we ask whether there exists a planar realization with one horizontal 
straight-line segment per color, called bus, such that all points with the same color 
are connected with vertical line segments to their bus. We present an ILP and an 
FPT algorithm for the general problem. For restricted versions of this problem, 
such as when the relative order of buses is predefined, or when a bus must be 
placed above all its points, we provide efficient algorithms. We show that another 
restricted version of the problem can be solved using 2-stack pushall sorting. On 
the negative side we prove the NP-completeness of a special case of BEP. 


1 Introduction 

Visualization of sets is an important topic in graph drawing and information visual¬ 
ization and the traditional approach relies on representing overlapping sets via Venn 
diagrams and Euler diagrams [32]. When more than a handful sets are present, how¬ 
ever, such diagrams become difficult to interpret and alternative approaches, such as 
compact rectangular Euler diagrams are needed [31], 

Often the geometric position of the elements of the sets are prescribed as points in 
the plane. The task is to emphasize the sets where the elements belong to. In visual¬ 
ization approaches for set memberships of items on maps, this is done by connecting 
points from the same set by corresponding lines (LineSets [2]), tree structures (KelpEu- 
sion [27]), and enclosing polygons (BubbleSet [11] or MapSets [13]). 

We consider a unified version of the tree-structure approach using a model that has 
been applied before for drawing orthogonal buses known from VLSI design [25, 34], 
Our goal is a membership visualization of points in sets by a tree-structure that consists 
of a single horizontal segment, called bus, to which all the points from the same set are 
connected by vertical segments, called connections', see Eig. 1 for planar and non-planar 
versions. We assume the sets to be given by single-colored points, such that in the final 
visualization, called bus realization, every point of the same color is connected to ex¬ 
actly one bus associated with this color. The objective is to find a position for each bus, 
such that crossings of buses with connections are avoided, called planar bus realization. 
We call this the bus embeddability problem (BEP). Such a simple visualization scheme 
makes it very easy to recognize the sets and label them, by placing a label inside each 
bus (if the bus is drawn thick enough), or directly above/next to the bus. 





(a) 


(b) 


(c) 


(d) 


Fig. 1. (a) Fixed positions of points, where points with the same color belong to the same set. 
(b) A planar bus realization for this setting, while (c) is a non-planar bus realization, (d) A point 
set without any planar bus realization. 


Related Work. Buses have been used, in a more general form, for visualizing degree- 
restricted hypergraphs. Ada et al. [1] used horizontal and vertical buses in bus real¬ 
izations, where the points (representing hypervertices contained in at most four hyper¬ 
edges) were not predefined in the plane. They asked whether a given hypergraph admits 
a non-planar bus realizations (allowing connections to cross each other) and showed 
that the problem is NP-complete. In contrast, if a planar embedding is given, a planar 
bus realization can be constructed on a 0{n) x 0{n) grid in time [7]. These 

types of problems also have connections to rectangular drawings, rectangular duals and 
visibility graphs, since the edges of the incidence graph of a hypergraph enforce visi¬ 
bility constraints in the bus realizations [19,33]. 

Another related approach is visualization based on graph supports of hypergraphs. 
Here the goal is to connect the vertices in such a way that each hyperedge induces 
a connected subgraph [6,8,23]. Supported hypergraph visualizations inspired edge¬ 
bundling and confluent layouts as alternative visualizations for cliques [12,15,30]. 

A solution to the BEP problem can be viewed as planar tree support for hyper¬ 
graphs, and this problem is related to Steiner trees [21], where the goal is to connect 
a set of points in the plane while minimizing the sum of edge lengths in the resulting 
tree; this is a classic NP-complete problem [16]. Hurtado et al. [20] considered planar 
supports for hypergraphs with two hyperedges such that the induced subgraph for ev¬ 
ery hyperedge and the intersection is a Steiner tree. Their objective was to minimize 
the sum of edge lengths, while allowing degree one or two for the hypervertices. BEP 
is even more closely related to rectilinear Steiner trees [14], where the Euclidean dis¬ 
tance is replaced by the rectilinear distance; constructing rectilinear Steiner trees is also 
NP-complete [17]. A single trunk Steiner tree [10] is a path which contains all vertices 
of degree greater than one. This is a variant that is solvable in linear time. BEP for 
a single set is the single trunk rectilinear Steiner tree problem, where we ignore the 
minimization of the sum of the edge lengths. Thus BEP can be seen as a simultaneous 
single-trunk rectilinear Steiner tree problem. The fact that a bus placement influences 
the placement of other buses makes the problem hard. 

Consider the input to BEP along with a box that encloses all the points. If in BEP 
the buses extend to the right boundary of this box, or both to the left and right bound- 
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ary of this box, then this problem corresponds to backbone boundary labeling and can 
be efficiently solved [4]. In backbone boundary labeling, the problem is to orthogo¬ 
nally connect points by a horizontal backbone segment leading to a label placed at the 
boundary. In this setting it is always possible to split the problem into two independent 
subproblems, which is impossible in our case. 

BEP is also related to the classical point set embeddability problem, where given 
a set of points along with a planar graph, we need to determine whether there exists 
a mapping of vertices to points such that the resulting straight-line drawing is planar. 
The general decision problem is NP-hard [9]. In the variant of orthogeodesic point set 
embedding, Katz et al. proved that deciding whether a planar graph can be embedded 
using only orthogonal edge routing is NP-hard [22]. 

Our Results. In Section 2 we solve BEP when the relative order of the buses is pre¬ 
scribed; we also show that BEP is fixed-parameter tractable (EPT) with respect to the 
number of colors. In Section 3 we formulate an integer linear programming (ILP) for¬ 
mulation for BEP and show some experimental results. In Section 4 we restrict BEP 
(when a bus must be above all its points, or a bus must be either at its topmost or bot¬ 
tommost point) and describe efficient algorithms for these settings. Another restricted 
version of the problem is shown to be equivalent to the problem of sorting a permuta¬ 
tion, which is called 2-stack pushall sorting. Einally we prove that BEP is NP -complete, 
even for just two points per color, if points may not lie on buses. 


2 Preliminaries 

We begin with some definitions. Suppose we are given a set of points V = {pi,... ,p„} 
and colors C = {ci, ..., Cfe} together with a function / : V — > C, /(p) = c. Eor simplic¬ 
ity, we assume that no two points share a coordinate in the input point set, although in 
some illustrations the input points might violate this assumption. The bus embeddabil¬ 
ity problem (BEP) asks, whether there is a planar bus realization with one horizontal 
bus per color. BEP is a decision problem, but in our descriptions whenever the answer 
is affirmative we also compute a drawing. We refer to such a drawing as a solution of 
BEP. In the negative case, we say that BEP has no solution. 

A point p has x-coordinate x{p), y-coordinate y{p), and color /(p). In a bus re¬ 
alization we have connections only between a point p and a bus c of the same color, 
that is, c = /(p). We denote by the set of points with color c. Bus c naturally 

extends from the x-coordinate xi{c) = min{x(p)|p e /"^(c)} of the leftmost point 
to the x-coordinate Xr(c) = max{x(p)|p 6 /“^(c)} of the rightmost point of /"^(c). 
We call 'ixi{c),Xr{c)] the span of c, which is predefined by the input points. The y- 
coordinate of a bus c is denoted by p(c), which is the only parameter to be determined 
for a solution for BEP. 

Note that BEP is trivial when there are at most two colors: it is always possible to 
place one bus at the top and the other (if exists) at the bottom of the drawing. Thus in 
the following we assume k> 2. Eor more than two colors, the relative order of the buses 
is important; see Eig. 1. Suppose the y-order of the buses is prescribed. The next lemma 
shows that one can check an existence of a solution for BEP respecting the order. 
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Lemma 1. There is a O{n\ogn)-time algorithm that, given an order of buses, tests 
whether there exists a solution for BEP respecting the order. 

Proof Suppose we are given an order ci < • • • < Cfe of the buses from bottom to top. We 
use discrete values for the y-coordinates increasing from bottom to top, where a unit is 
1/n of the y-distance of two consecutive points. We first present a simpler 
algorithm, and then describe how to speed it up. 

Recall that the span of every bus is defined by an input point set; hence, we only 
show how to choose y-coordinates of the buses. The first bus, ci, is placed at y-coordinate 
y{ci) = 0, and all the points of color ci are connected to the bus. Assume that bus 
Ci_i is placed at y-coordinate y{ci-i) and is connected to all its points. We place 
at y{ci) = y{ci-i) + 1 unit and check if the bus crosses a previously drawn (vertical) 
segment. If it does cross a segment, then we shift Ci one unit upwards by increasing 
y{ci) and repeat the procedure. Once the bus is placed without crossings, we connect it 
to the corresponding points. Consider the vertical segment of a point p of color Ci. It is 
easy to see that if y{p) > y{ci), then the segment cannot cross a previously placed bus 
Cj for j < i. If y{p) < y{ci) and the vertical segment crosses a bus, then such a cross¬ 
ing is unavoidable in any solution respecting the given order. Hence, we may stop the 
algorithm reporting that no solution exists. Otherwise, we proceed with the next color. 

The above algorithm can easily be implemented in quadratic time. However, we can 
do better using the following observation: Every bus is placed at its bottommost “valid” 
y-coordinate, that is, the one that does not produce crossings with previously placed 
buses. To find such a y-coordinate efficiently for each color, we store all points of the 
already processed colors in a data structure D that supports the range operation such 
as “extracting minimum/maximum on a given range”. For every color q, we extract a 
point with the maximum y-coordinate in the range corresponding to the span of q. The 
bus of Ci is placed at the maximum of the extracted y-coordinate and the y-coordinate 
of bus y{ci-i). Then all the points of color q are added to D. A balanced tree (e.g., 
a segment tree) providing logarithmic complexity for insert and extract operations is 
sufficient for our needs. □ 

In general the correct order of the buses for a planar bus realization is not known. 
One can apply Lemma 1 for each of the k\ possible bus orders, which yields an (D{k\)- 
time^ algorithm for BER Next, we improve the running time with an algorithm provid¬ 
ing deeper insight into the structure of the problem. 

Lemma 2. There is a 0{2^)-time algorithm for BEP. 

Proof. We solve a given instance of BEP using dynamic programming. Let us call 
a state a pair {h,B), where 0</i<n-i-lisan integer and B is a subset of C = 
{ci,..., Cfe}. By a solution for a state {h,B) we mean a (planar) bus realization consist¬ 
ing of buses for every color ce B such that the topmost bus has y-coordinate h. If such 
a solution exists, we write F{h, B) = true, and otherwise F{h, B) = false. It is easy to 
see that a solution for the original BEP problem exists if and only if F{h,C) = true for 
some 0 < ft. < n H-1. 

* O hides polynomial factors. 
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We reduce the problem to solving it for “smaller” states, that are the states with 
fewer elements in B. As a base case, we set F{h,B) = true for all 0 < ft. < n + 1 
and \B\ = 1. To compute a value for a state F{h, B) with \B\ > 1, we consider a color 
c* 6 B. Let ft* = max{y(p)|/(p) e B \ {c*} and xi(c*) < x{p) < Xr{c*)}, that is, the 
largest (topmost) y-coordinate of a point of color B \ {c*} laying in the span of c*. It 
follows from the proof of Lemma 1 that the bus for c* should be placed at y-coordinate 
ft*. Thus, F{h, B) is set to true if (a) ft > ft* and (b) there exists a solution for a state 
(ft', B \ {c*}) for some ft' < ft. We stress here that in order to compute F{h, B), one 
needs to consider every color of S as a potential c*. There are n2^ different states, and 
a computation for a single state clearly takes a polynomial number of steps. □ 

The above result shows that the BEP problem is fixed-parameter tractable with re¬ 
spect to k, that is, it can be efficiently solved for a small number of buses. Note that in 
Section 5 we prove that BEP is NP-complete; hence, it is unlikely that a polynomial¬ 
time (in terms of k) algorithm exists. 

3 An ILP for BEP 

In this section we present an integer linear programming (ILP) formulation for BEP that 
produces a planar bus realization if one exists. The ILP also minimizes the amount of 
ink in a solution, that is, the sum of all segment lengths. 

Lemma 3. A solution for BEP can be computed by an ILP. 

Proof. In a preprocessing step we compute the span of every bus c e C. As mentioned 
earlier, it remains to compute the y-coordinate variable y{c) of every bus c. To this end, 
we introduce a planarity constraint for every point p eV within the span of bus c having 
a different color. The pairs {p,c),c+ f{p) are called conflicting. Conflicting pairs {p, c) 
are stored in a matrix J and induce the constraint {y{p) < y{c) and y{f{p)) < y{c)) 
or {y{p) > y{c) and y{f{p)) > y{c)). The matrix J' can be computed in 0{kn) time, 
where n=\V\ and k=\C\. In order to minimize the amount of ink, we sum up the lengths 
of all connections and ignore the lengths of buses, as those are determined by the input. 

min 

“C f{p)=c 

s.t. {y{p)<y{c)\/y{f{p))>y{c))A{y{p)>y{c)\/y{f{p))<y{c)) 'i{p,c)ej 
0 < y{c) < max{y(p)} + 1 

p^V 

Since absolute value (resp. “or”) needs one more variable and 3 constraints for ev¬ 
ery point (resp. for every conflicting pair) the final ILP has n+k+2\J'\ variables and 
3n+k+6\J^\ constraints. □ 

In order to get a feeling about the probability that a point set admits a solution of 
BEP, we ran an experiment with the ILP, implemented with the Gurobi solver [18]. We 

^ min X! |a - i>| min f, e, s > a - b, e > b - a, e > 0; {a < b) v (c < d) <=> a - b < eM, c- d < 
(1 - e)M, e 6 {0,1}, M = oo 
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Distribution of solutions for BEP 
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Fig. 2. The percentage of solutions for BEP for a random point set of size n = kl with I = 2,3,4 
points per color out of fc = 3,..., 20 colors. 

considered point sets with k = 3,... ,20 colors and with Z = 2,3,4points per color. We 
randomly placed the points on a 1024 x 768 area. For each pair (/, k) we counted the 
number of BEP solutions out of 100 instances; see Fig. 2. The remaining instances were 
infeasible. For a fixed number of points, I, the number of solutions for BEP decreases 
with increasing the number of colors, k. It decreases faster the higher I is. On the other 
hand for a hxed number of colors, k, the number of solutions for BEP also decreases 
with increasing number of points, 1. Hence, studying two points per color promises to 
be sufficiently interesting. Thus, as the base case for further analysis, we initially con¬ 
sider two points per color, before dealing with the general case, where in real instances 
solutions rarely exist. It is possible that much more solutions exist if we allow only few 
crossings, but all non-planar settings are left as open problems. 

4 Efficiently Solvable BEP Variants 

In this section we consider three variants of BEP, which can be solved in polynomial 
time. A bus c is called top (resp., bottom) if all of its points are below (resp., above) 
the bus, that is, y{c) > y{p) (resp., y{c) < y{p)) for all p 6 f~^{c). We distinguish 
between buses that are above (below) of their points and buses that pass through one of 
their points. A top-bus is a n-bus if y{c) > y{p) for all p 6 /“^(c) (Fig. 3(a)), while it 
is a r -bus if y{c) = y{p) for a point p with y{p) = max{y{q)\q € /“^(c)} (Fig. 3(c)). 
Similarly we dehne a u-bus and a \--bus-, see Figs. 3(b) and 3(d). A bus, whose type 
is none of the four types from above, is called a center-bus. The variant of BEP where 
only buses of the types in S' £ {n, u, r, l} are allowed to use is denoted by S-BEP. 

In Section 4.1 we study n-buses and provide an algorithm for n-BEP. The same 
algorithm obviously solves the u-BEP variant. Next we consider r-buses and L-buses. 
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Note that r-BEP and l-BEP are trivial, since every r-bus (resp., L-bus) is uniquely 
defined by its span and the topmost (bottommost) point. Hence, we investigate and 
design an efficient algorithm for the (r,L)-BEP variant. Einally in Section 4.3, we 
examine the general BEP for a specific point set, where all points lie on a diagonal. 
We show that the variant of the problem is equivalent to a longstanding open problem 
(resolved very recently) of sorting a permutation with a series of two stacks. 

4.1 n-BEP 

Here, we present an algorithm that decides in polynomial time whether a drawing with 
n-buses exists for a given input, and constructs such a drawing if one exists. 

Theorem 1. There exists an 0{n\ogn)-time algorithm for n-BEP. 

Proof. Eor ease of presentation, we first assume that the input consists of two points 
per color, that is, k = n/2, and provide a simple quadratic-time implementation. Later 
we generalize the algorithm and improve the running time. Intuitively, the algorithm 
sweeps a line from bottom to top and processes the points in increasing order of y- 
coordinates. At every step, we keep all the vertical segments of the “active” colors (the 
ones without a bus) in the correct left-to-right order. If two vertical segments of the same 
color are adjacent in the order, then we can draw the corresponding bus and remove the 
color and its vertical segments. Otherwise, all the active vertical segments have to be 
“grown” until we reach the next point. It is easy to see that a solution exists if and only 
if the set of active colors is empty after processing all the points. 

More formally, the points are processed one-by-one in increasing order of their y- 
coordinates. The points are stored in an array sorted by x-coordinate, that is, we have 
{pi,... ,pn) with x{pi) < ••• < x{pn). At each iteration, a new point is inserted into 
the array in the position determined by its x-coordinate. Then the array is modified 
(or simplified) so that the pairs of points of the same color that are adjacent in the 
array are removed. That is, if f{pi) = f{pi+i) for some 1 < i < n, then we get a 
new array (pi,... ,Pi-i,Pi+ 2 , ■ ■ ■ ,Pn)- The simplification is performed as long as the 
array contains monochromatic adjacent points. After this step the algorithm proceeds 
with the next point. Eor every color c, we keep the value y*{c), which is equal to the 
y-coordinate y{p),p € f~^{c') of the point of color c', whose insertion into the array 
induced the removal of points /"^(c) from the array. If the algorithm ends up with a 
non-empty array, then we report that no solution exists. Otherwise, the y-coordinate 
of the resulting bus of color c is y*{c) + e, where e > 0 is sufficiently small to avoid 
overlaps between the buses. An example of the algorithm is illustrated in Eig. 4. 


(a) (b) (c) (d) 

Fig. 3. Illustration of (a) n-bus, (b) u-bus, (c) r-bus, and (d) L-bus. 
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Steps: 


Array: 

8 RWRW 
7 RBBWR —RWR 

6 RBWR 
5 RBGGR — RBR 

4 RBGR 
3 RGR 
2 RG 
1 R 


Fig. 4. Running the algorithm from Lemma 1 on a given point set with red (R), green (G), blue 
(B), and white (W) pairs of points. Since the resulting array is not empty, there is no solution for 
the instance. Notice that removing any of the colors yields an instance with a solution. 


Correctness. The correctness follows from the observation that the algorithm chooses 
the lowest “available” y-coordinate for every bus, that is, the one that does not induce 
a crossing between the bus and vertical segments of other colors. Indeed, if at any step 
of the algorithm we get a color pattern R,..., B,..., Rin the array formed by red (R) 
and blue (B) points and the second blue point p has not been processed yet, then clearly 
in any solution the red vertical segments reach the y-coordinate of p. Hence, it is safe to 
“grow” the segments. On the other hand, if processed points form a color pattern RR 
(that is, two consecutive points of the same color), then there is a solution connecting 
the corresponding vertical segments at the current y-coordinate. The two points can be 
removed from consideration, as they cannot create crossings with the subsequent buses. 
It is also easy to see that the algorithm minimizes ink of the resulting drawing. 

Running time. At every iteration of the algorithm, we need to insert a new point into the 
sorted array and then run the simplification procedure. Point insertion takes 0{n) time 
and removal of a pair of points from the array can also be done in 0{n) time. Since 
every pair is removed only once, the total running time is 0{n^). 

To get down to O{n\ogn) time, we use a balanced binary tree instead of an array 
to store the points. The tree is sorted by the x-coordinates of the points; hence, inser¬ 
tion/removal of a point takes O{\ogn) time. Note that after inserting/removing a point, 
the only potential candidate pairs for simplification are the point’s neighbors that can 
be found in O{\ogn) time. Again, every point is inserted/removed only once; thus, the 
total running time is 0{n log n). 

Finally, we observe that the algorithm can be generalized to handle multiple points 
per color. To this end, we change the simplification step so that the points are removed 
only if they form a contiguous subsequence in the array (tree), containing all points of 
this color. Hence we need to know the number of points for each color, which can be 
done with a linear-time scan of the input. It is easy to see that the proof of correctness 
can be appropriately modified and the running time remains the same. □ 

4.2 (r,L)-BEP 

We present an algorithm that decides in polynomial time whether (r,L)-BEP has a 
solution for a given input, and constructs a drawing if one exists. 

Theorem 2. There exists an 0{n^)-time algorithm for (r, \_)-BER 
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X, = true 




Xq => X, 


Fig. 5. Three examples for creating clauses for two colors black and white. 


Proof. The span of every bus is predefined by the input, while the y-coordinate has 
precisely two options. We show that (r, l)-BEP can be modeled by 2-SAT, and thus is 
efficiently solvable. For ease of presentation, we first assume that the input consists of 
two points per color and describe a simple quadratic-time algorithm. 

The algorithm creates a variable Xc for every color ceC. The value of Xc is true if c 
is a r -bus, and it is false if c is a L-bus. Then for every pair of colors c, c', the algorithm 
creates a clause for the 2-SAT instance when the corresponding buses induce a crossing. 
Building the clauses with respect to the relative position of points is a straight-forward 
procedure; 3 examples are illustrated in Fig. 5. 

Specifically we create clauses according to the following analysis. Fet R{c) be 
the smallest enclosing rectangle of the points pc, Qc of color c. By symmetry, we may 
assume that pc appears in the left bottom corner, while q,. appears in the right top corner 
of R{c). 

We distinguish the cases when 

(1) points Pc',qc' are in the top left, bottom right corner of R{c') or whether 

(2) points Pc',qc' are in the bottom left, top right corner of R{c'). 

In each of the two cases we consider the 8 subcases, which are 

(a) R{c') intersects only the top boundary of i?(c), 

(b) R{c') intersects only the bottom boundary of R{c), 

(c) R{c') intersects only the right boundary of i?(c), 

(d) R{c') intersects only the left boundary of i?(c), 

(e) R{c') contains the top right corner of i?(c), 

(f) R{c') contains the bottom right corner of i?(c), 

(g) R{c') contains the top left corner of i?(c), 

(h) R{c') contains the bottom left corner of R{c). 


cases 

a 

b 

c 

d 

e 

f 

9 

h 

1 

Xc = f 

Xc=t 

Xc=t 

Xc = f 

Xc' = t 

Xc = t 

Xc = f 

Xc' = f 

2 

Xc = f 

Xc=t 

Xc=t 

Xa = f 

Xc'^Xc' 

Xc' => Xc 


Xc^Xc' 


Table 1. For each of the cases (a)-(h) from above we build a clause depending on the configuration 
(l)-(2) from above, where t stands for true and / for false. 


Correctness. The correctness follows from the complete case analysis by the rules of 
Table 1. 
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Running time. We remark that for the n^/4 pairs of colors, we create (!?(n^) clauses, 
each clause in constant time by a case analysis. This results in a 2-SAT instance with 
k variables Xc, c e C and 0{n^) clauses. We solve this instance in linear time [3] and 
the solution determines the drawing: c is drawn as a r-bus, if the value of Xc is true, 
otherwise c is drawn as a L-bus. 

We can generalize this idea to the case of more points per color. In the general case 
the y-coordinate of a bus again has precisely two options. In contrast to the case with 
two points per color we check several points (not only the leftmost or rightmost point) 
of color c' for their position with respect to the points of color c, since points lie not 
necessarily in corners of the enclosing rectangle. □ 


4.3 Diagonal BEP 

Here we consider a diagonal point set in which all points lie on a single diagonal line 
and there are two points per color. We assume that the point set is separable, that is, 
there is a straight line separating every pair of points having the same color; see Fig. 6 . 
This specific arrangement can be naturally described in terms of permutations. Assum¬ 
ing that the colors are numbered from 1 to /c in the order along the diagonal from bottom 
to top, the input is described by a permutation tt = [ 7 r(l),..., 7 r(A:)] on {1,..., fc}. Such 
an instance is called diagonal tt-BEP. 

It turns out that this variant of BEP is closely related to the well-studied topic of 
sorting a permutation with stacks introduced by Knuth in the 1960’s [24]. We next 
show that diagonal tt-BEP has a solution if and only if tt can be sorted with 2 stacks in 
series. The problem of deciding whether a permutation is sortable with 2 stacks in series 
is a longstanding open problem and it has been conjectured to be NP-complete several 
times [5]. Only very recently a polynomial-time algorithm has been developed [28,29]. 
It is an indication that even the restricted variant of BEP is highly non-trivial. Next we 
prove the equivalence. 

Eirst observe that for a diagonal point set with 2 points per color, a top-bus (bottom- 
bus) can be transformed to a center-bus. Eor every color c, there are no points of different 
color within the span of c above the topmost point of c. Hence, we may only consider 
center-buses in the variant of BEP. 

Eor the 2-stack sorting problem, given a permutation tt, we want to sort the num¬ 
bers to the identity permutation [1,..., fc] with two stacks Sj, Su using the following 
operations: 

• ai : read the next element i from input tt and push it on the first stack Sj', 

• fii : pop the topmost element i from Si and push it on Sif, 

• : pop the topmost element i from Su and print it to the output. 

To proof of the equivalence between 2-stack sorting and bus embeddability, we note 
that the first operation, ai, corresponds to the left vertical segment of color i, the second 
one, Pi, is the bus of i, while 7 ^ corresponds to the right vertical segment of the color; 
see Eig. 6 and Eig. 7. A crossing in the drawing correspond to an “invalid” sorting 
operation in which either a non-topmost element is moved from Si to Su (a crossing 
to the “left” of the diagonal), or a non-topmost element is moved from Su to the output 
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Fig. 6. A diagonal point set with a solution for BEP and the regarding sorting sequence. 


(a crossing to the “right” of the diagonal). Hence, sorting sequences of the operations 
for TT are in one-to-one correspondence with planar bus realization for the point set. 
Since the point set is separable, all the elements of tt will be pushed to Sj before any of 
the elements is popped to the output. This is called 2-stack pushall sorting [28] and is 
considered next in more detail. 

We can describe a sequence of the operations by a word w € {a, ^, 7 }^", where 
every operation appears n times. 

For example w = Q; 3 a 2 cricr 4 / 34 / 3 i 7 i/ 3272/337374 is a sorting word for tti = 3214 to 
712 = 1234 with two stacks, see Table 2. 


operation 

input 

Si 

Sii 

output 


3214 




Q3 

214 

3 



OL2 

14 

23 



Ql 

4 

123 



0L4 


4123 



Pa 


123 

4 


Pi 


23 

14 


71 


23 

4 

1 

p2 


3 

24 

1 

72 


3 

4 

12 




34 

12 

73 



4 

123 

74 




1234 


Table 2. Permutation [3,2,1,4] is sortable with two stacks. 


A word w also encodes the input and output of a sequence by subscripts, when 
disregarding the subscripts of the beta operation. For example s{w) = 32141234 is the 
sequence of subscripts for w. 
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output 


Fig. 7. A correspondence between 2-stack sorting and a planar bus realization. 


We may restrict this sequence of operations to only a and 7 operations, denoted 
by tt;|{Q;, 7 }. We say w is a pushall word, if s(r(;|{a, 7 }) = 7ri7r2. The word w' = 
a3a2ai/3i7iQ;4/34/?272/?37374 also sorts tti to 772 with two stacks, but w' is not a 
pushall word, since s (w') = 32114234 7 7ri7r2. 

Now we assume we are given 2n points on a diagonal respecting the order tti,tt 2 . 
We denote by 7ri(7r2) the order of the first (second) appearance of the elements. Every 
output word w of the 2-stack-pushall-sortable algorithm describes the sorting from tti 

to 742. 

The 2-stack-sorting algorithm takes as input tti and 742 and returns in 0{ri?) time 
one sorting word of E = {w ■ w sorts tti to 742 }. If such a word w exists, then we 
can construct a planar bus realization with center buses of the embedded points 7ri7r2 
according to w as follows. We apply one of the following 3 rules on the letters of w. We 
process w letter by letter and read along 3n imaginary slots on the diagonal. 

ai the next slot of the diagonal is point i with a connection going up. 

Pi the next slot of the diagonal is taken by the horizontal segment from the end of the 
connection of point i, then crossing the diagonal. 

7i the next slot of the diagonal is point i with a connection down to its horizontal 
segment, extended such that this connection meets perpendicular. 

This drawing is planar: 

• any crossing of two edges incident to i,j to the left of the diagonal comes from the 
sequence ..., a(i), ..., a{j), ..., /3(i), ..., which means push i on Sj, then push 
j on Si and then pop i from Si, which is impossible since i is not the topmost 
element of S'/. 

• any crossing of two edges incident to i, j to the right of the diagonal comes from the 
sequence ..., /3(i), ..., P{j), • • ■, j{i), ■ • •, which means push i on S//, then push 
j on Sii and then pop i from S// (and print i to the output), which is impossible 
since i is not the topmost element of S//. 

The construction from a planar bus realization with center buses of a diagonal point 
7ri7r2 set to a sorting word w for 744742 is just traversing the diagonal from bottom to top 
and simultaneously building incrementally the sorting word w. We start with w = X, 
where A is the empty word. If the next item on the diagonal is the first appearance of 
a letter i, we set w = ru o ai. If the next item on the diagonal is a crossing of the edge 
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Fig. 8. A clause, variable and chain gadget for reduction from planar 3-SAT. Vertical propagation 
of true and false are unique, but n-buses are just uniquely propagated in the top direction and 
u-buses are just uniquely propagated in the bottom direction. 


connecting the two points of i, we set tt; = w o . If the next item on the diagonal is the 
second appearance of a letter i, we set re = re o . It is easy to see that this word w sorts 
TTi to tt 2 - This finishes the proof. 

Theorem 3. Diagonal tt-BEP has a solution if and only if n is 2-stack pushall sortable. 
This can be checked in O(n^) time. 

5 Hardness of BEP 

In this section we consider BEP^, where e > 0 is the additional input number indicating 
the minimum allowed distance between points and their bus. We prove that BEP® is 
NP-complete even for 2 points per color. 

We can easily verify a possible solution using Lemma 1; thus BEP® is in the class 
NP. We then show that (n,u)-BEP® for 2 points per color is NP-hard. To prove the 
hardness of (n,u)-BEP®, we reduce from planar 3-SAT [26], which is 3-SAT, where an 
instance is represented by a graph whose vertices represent variables and clauses and 
whose edges represent containment of variables in clauses. The most important module 
of the construction is a chain link, which is also a gadget for replacing variables. It 
consists of two points on a common horizontal line that will be connected by a bus. We 
replace the edges of the graph by chains consisting of nested chain links and replace 
the clause vertices by a big construction of points, that allows two specific points to 
be connected via a bus using only one of three choices, cf. Pig. 8. We use the input 
£ to be able to block some choices for this bus. We first restrict ourselves to (n, u)- 
BEP® and drop the “no points share a coordinate” restriction. We finally transform the 
construction into the “no points share a coordinate” setting and allow also center-buses. 

A variable gadget consists of two points oi, 02 of the same color on the same y- 
coordinate. The value of the variable is true if the two points are connected with a n-bus 
and the value of the variable is false if the two points are connected with a u-bus. We 
use a variable gadget, referred to as a chain link, also as elements of chain gadgets. 
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A chain gadget propagates the value of a chain link, which is actually a variable 
gadget, to another chain link. Let Ui, 02 (respectively bi, b^) be the two points of the 
chain link at the beginning (respectively end) of the chain. A chain gadget consists of k 
chain links. 

In a horizontal chain gadget we place the points on a single horizontal line in the 
order oi, 5i, 02 , &2 (respectively 5i, oi, & 2 , ^ 2 ) for propagating to the right (respectively 
to the left). If oi, 02 are connected with a n-bus, then bi, 62 must be connected with a u- 
bus and the other way round. This construction can be repeated until the chain consists 
of k chain links. The sign of a horizontal chain is defined by Clearly if the 

sign is positive, then the first bus and the last bus are of the same type. If the sign is 
negative, then the first bus and the last bus are different. 

In a vertical chain gadget we place bi below (respectively above) oi and &2 below 
(respectively above) 02 on the same x-coordinate with a distance of 2e for propagating 
to the top (respectively to the bottom). It is easy to check that in such a way we can 
only uniquely propagate a n-bus to the top and a u-bus to the bottom. It may happen 
that the type of buses change during a vertical propagation. The sign of a vertical chain 
is defined as + 1 . 

The sign of two chains, which are connected, will be multiplied. If a literal in a 
clause appears positive, then the corresponding chain has sign - 1 , otherwise + 1 . 

A clause gadget consists of two main points pi,pr, 4 horizontal bounding segments 
SthStr, sujSbr, 8 Vertical bounding segments si ^,... , 841 , 54 ^, and 18 chain links, 
see a schematic illustration in Fig. 8. We aim at satisfying the clause if and only if a bus 
connecting the main points can be drawn. 

Within a bounding square Q we place horizontal bounding segment Sti (str) in the 
top left (right) corner, and horizontal bounding segment su (s&r) in the bottom left 
(right) corner. Above su (sbr) we place main point pi ipr), such that there is a normal 
to Sbi (sbr) through pi (pr) that is also crossing sti (str)- This construction prevents the 
bus connecting the main points to be in the exterior of Q. 

In Q there are two vertical lines I and r that both separate pi from pr. We place the 
vertical bounding segments six,S 2 x,S 3 x,S 4 x,x e {l,r} in this order from bottom to 
top on line x with e distance between every pair of consecutive segments. The resulting 
horizontal space between segment Six and Si+xx is called i-th gap, i = 1,2,3. The gaps 
represent the literals in the clause. This construction restricts the choices for the bus 
connecting the main points to be precisely three. 

Finally we place 9 chain links below the first gap, 6 chain links between the first and 
second gap and 3 chain links between the second and the third gap. More specifically 
let ui,..., ue be 6 vertical lines between I and r in this order from left to right. We place 
3 chain links on lines vi , V 2 such that the first chain link has its points on the boundary 
of Q, the last chain link has its points on the bottom boundary of the first gap and the 
distance between every pair of chain link points is at most 2£. Similarly we place 6 
chain links on lines V 3 , V 4 such that the first chain link has its points on the boundary of 
Q, the 4th chain link has its points on the top boundary of the first gap, the last chain 
link has its points on the bottom boundary of the second gap, and the distance between 
every pair of chain link points is at most 2e. In the same way we place 9 chain links 
on lines U 5 , ue such that the first chain link has its points on the boundary of Q, the 4th 
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(7th) chain link has its points on the top boundary of the first (second) gap, the last chain 
link has its points on the bottom boundary of the third gap, and the distance between 
every pair of chain link points is at most 2e. We refer to the last chain link of lines 
vi , V 2 (respectively lines vz,Vi and V 5 ,vq) as the chain link of the first gap (respectively 
second and third gap). This construction allows to block or open gaps from the bottom 
of Q. 

Notice that it is easy to simulate vertical or horizontal segments with points as 
demonstrated in Fig. 9. 

The construction of an instance of (n,u)-BEP® from an instance / of planar 3- 
SAT is done according to a planar drawing of the graph Gj. We may assume that all 
variable vertices of Gi are on a single horizontal line. We use this line and place the 
variable gadgets according to this order. The clause vertices above the variable vertices 
(top clauses) are replaced by clause gadgets and the clause vertices below the variable 
vertices (bottom clauses) are replaced by horizontally mirrored clause gadgets. Finally 
we replace the edges of Gi with chain gadgets. 

The number of points needed to construct an instance of (n, u)-BEP® is polynomial 
in n and m. Given a planar 3-SAT instance with n variables and m clauses, the cor¬ 
responding (n,u)-BEP® instance has at most 0{nm) points. Eor a clause gadget we 
need precisely 118 points, for a variable gadget we need 2 points and for chain gadgets 
we need 10 points plus the points needed to surround other clause gadgets. If an edge 
from a clause to variables vertically passes k other clauses, then we need 18A: points 
for this construction. Since we have 0{n + m) edges and m clauses, we might need 
0 {nm + rrf) points for the edges. 

Pig. 9 shows how to use a variable several times: we stretch one chain link for 
horizontal propagation and add a vertical chain gadget for vertical propagation. 

Theorem 4. (n, u)-BEP^ for 2 points per color is -complete. 

Proof To show the membership of (n,u)-BEP^ in the class NP we observe that we 
have n points and between every pair of consecutive points we have a gap. In every gap 
there can be possibly n buses, that is, we have (n - l)n slots, where to place buses. So 
every slot represents a possibility to place a bus. We can guess a drawing by choosing an 
order of the buses: all the drawings where buses move within their gap are equivalent. 
To check if the order leads to a feasible solution of (n, lj)-BEP^, we apply the algorithm 
of Lemma 1. 

We prove the hardness of (n, lj)-BEP^ by a reduction from planar 3-SAT [26]. Let I 
be an instance of the planar 3-SAT problem and let Pj be the point set constructed from 
the gadgets, that is, we replace in the planar graph representing I every clause vertex 
by a clause gadget, every variable vertex by a variable gadget and every edge by a chain 
gadget. We prove next that Pj admits a solution of (n, lj)-BEP^ / has a satisfying 
truth assignment. 

If Pj admits a solution, then in particular every pair of main points is con¬ 
nected. Consider w.l.o.g. a top clause c with literal y corresponding to the gap through 
which the main points are connected. We associate with y the gap of c. If the chain link 
of y is a u-bus, then this bus is uniquely propagated to the bottom and does not change 
its type. If y is a positive variable x, then by construction the chain has sign -1 and 
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Fig. 9. Point set instance constructed via gadgets for I = {xi v X2 v xs) A {x2 v X3 v X4). The 
red buses indicate the truth assignment X3 = X4 = f, Xi = X 2 = t. Clause gadgets are enclosed by 
a green rectangle. The green thin circles indicate that distances are less than 2e, while the green 
thick circle indicates a change of bus types during a not-unique vertical propagation of top buses 
to the bottom, that is, X 3 is meaningless for the second clause. 


the chain ends in a n-bus at the variable gadget, corresponding to x being true. If y is 
a negated variable x, then by construction the chain has sign +1 and the chain ends in 
a u-bus at the variable gadget, corresponding to x being false. For bottom clauses the 
same argument holds horizontally mirrored. 

Assume we have a satisfying truth assignment for I. We explain how to con¬ 
struct a solution for (n,u)-BEP'^. First for every variable being true we draw a n-bus, 
while for every variable being false we draw a u-bus. We propagate n-buses with n- 
buses to the top and to the bottom, while we propagate u-buses with u-buses to the top 
and to the bottom. A n-bus (u-bus) ends in a n-bus (u-bus) if the variable in the top 
clause is negated (positive) or the variable in the bottom clause is positive (negated). 

We keep the type of buses in a vertical propagation as long as possible, which can 
only be interrupted by a main bus. Then we change the type of buses and the gap be¬ 
comes blocked, although the variable is true and appears positive, or the variable is false 
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and appears negative in the clause. Additionally this interrupting main bus indicates that 
this clause is already satisfied and thus the variable of the interrupted chain is irrelevant 
for the satisfiability of this particular clause. 

By construction we get a feasible (planar) solution for the buses in P/. 

Finally we translate the construction into the “no points share a coordinate” set¬ 
ting. We may assume an underlying k x k grid with grid unit e/2 in the plane so 
that all points have integer coordinates. Let pi,..., be the points ordered first by x- 
coordinate, then by y-coordinate. We modify the x-coordinates by a shift x{l) = a;(() + l 
for all I > j, if x{pi) = x(pj), as long as two points share the same x-coordinate. We 
apply the same modification in y-direction with respect to the same order of points 
Pi,. ■. ,Pn- Finally no points share a coordinate. 

The properties of depending colors stay the same, since the topological operation is 
just a stretch. Clearly chain links are dependent before the stretch, if and only if they 
are dependent after the stretch. □ 

We consider as an example the instance I = {xi VT2 v X3) a (x2 v Ta v a;4) of planar 
3-SAT. Clearly the clause-variable graph is planar. Fig. 9 illustrates the point set created 
from the instance I. 

We can adopt the same construction when additionally using center buses. Now 
some vertical segments can be modeled by using just two points. In a clause gadget, we 
move one of the main points from bottom to top such that the bus connecting the main 
points is necessarily a center bus. The remaining parts are the same. Also for center 
buses we need e as input for the minimum distance of buses to their points. Notice that 
a bus c and a point p of different color c + c{p) may be closer than e, as well as two 
buses c, c' may be closer than e. 

Theorem 5. BEP^ for 2 points per color is HP-complete. 


6 Conclusion and Future Work 

We studied bus embeddability, where a set of colored points is covered by a set of 
horizontal buses, one per color and without crossings. We described an ILP and an 
FPT algorithm for the general problem and presented polynomial-time algorithms for 
several restricted versions. The general problem is shown to be NP-complete even for 
two points per color when points may not lie on buses. 

It is still open to determine the complexity of BEP in the following cases: 

• BEP using only center-buses; 

• (n, u)-BEP, that is, BEP without center-buses; 

• diagonal BEP with more than 2 points per color; 

• general BEP (in our construction, we use an extra e as a parameter). 

A natural generalization would be to allow both horizontal and vertical buses, as 
in [1,7]. Another variant might be to consider multi-colored points, where a point has 
to be connected either to all the buses of its corresponding colors, or to at least one of 
them. Eor point sets that have no solution for BEP with only one bus per color, we may 
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allow more than one bus or bound the number of crossings. Possible objectives in these 
scenarios are to minimize the total number of buses over all colors, to minimize the 
total number of buses, or to minimize the total number of buses if each tree can connect 
< k unicolored points. These objectives are even interesting if a solution to BEP exists. 
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